home *** CD-ROM | disk | FTP | other *** search
/ Just Call Me Internet / Just Call Me Internet.iso / com / computer / casio_st / prog_fx / oisac2 / l_s6440.txt < prev    next >
Encoding:
Text File  |  1994-03-01  |  809 b   |  18 lines
  1. 10 GOSUB400:ANGLE1:IFp>1 ORp=<0 THEN*
  2. 15 IFp=1 THENx=0:GOTO80
  3. 20 IFm=1 THENm=p:p=.5-p/2:GOSUB"LIB0:S6430":GOSUB400:p=m:m=1:x=1/y^2:GOTO80
  4. 30 IFn=1 THENn=m:p=p/2:GOSUB"LIB0:S6430":GOSUB400:n=1:p=p*2:x=y^2:GOTO80
  5. 40 IFm=2 THENp=1-p:m=n:n=2:GOSUB100:p=1-p:n=m:m=2:x=1/x:GOTO80
  6. 50 IFn>m THENp=1-p:d=n:n=m:m=d:GOSUB200:x=1/x:d=m:m=n:n=d:p=1-p:GOTO80
  7. 60 GOSUB200
  8. 80 y=x:RETURN
  9. 100 GOSUB"LIB0:S6420":GOSUB400:x=y:a=n-2
  10. 110 x=x/n*(1+((x-a)/2+(((4*x-11*a)*x+a*(7*n-10))/24+(((2*x-10*a)*x+a*(17*n-26))*x-a*a*(9*n-6))/48/m)/m)/m):RETURN
  11. 200 GOSUB100
  12. 210 d=x:c=p:GOSUB"LIB0:S6240":GOSUB400:p=c
  13. 220 z=n+m:z=EXP((z*LN(z/(n*x+m))+(n-2)*LNx+LN(n*m/z)-LN(4*PI)-(1/n+1/m-1/z)/6)/2):x=x+(y-p)/z
  14. 230 IFABS(d-x)>3E-4 THEN210 ELSERETURN
  15. 400 ONERRORGOTO500:RETURN
  16. 500 IFERR=1 THENCLS:ONERRORGOTO0
  17. 510 sr=1:RESUME80
  18.